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Abstract 

In the framework of the fey-factorization approach, the isolated prompt photon pair pro- 
duction in pp and pp collisions at high energies is studied. The consideration is based on the 
quark-antiquark annihilation, quark-gluon scattering and gluon-gluon fusion subprocesses, 
where the non-zero transverse momenta of incoming partons are taken into account. The 
unintegrated quark and gluon densities in a proton are determined using the Kimber-Martin- 
Ryskin prescription. The numerical analysis covers the total and differential production cross 
sections and extends to specific angular correlations between the produced prompt photons. 
Theoretical uncertainties of our evaluations are studied and comparison with the NLO pQCD 
calculations is performed. The numerical predictions are compared with the recent experi- 
mental data taken by the D0, CDF, CMS and ATLAS collaborations at the Tevatron and 
LHC energies. 
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1 Introduction 



The production of prompt photon pair^J in hadronic collisions at high energies is of great 
interest at present p~H5]. It provides a major background to searches for rare or exotic pro- 
cesses in both the Tevatron and the LHC In particular, it represents a major source of large 
and irreducible background to searches for the Higgs boson in the low mass range, where its 
first indications (with mass of about 125-126 GeV) have been found very recently. It is also 
a significant background to searches for a new heavy resonances [6] , cascade decays of a new 
heavy particles [7] and to searches for some effects of physics beyond the Standard Model. 
For example, so called Universal Extra-Dimensions predict non-resonant diphoton produc- 
tion associated with significant missing transverse energy [8j|9]. Other extra-dimension mod- 
els, such as Randall- Sundrum model [ID] , predict the production of gravitons, which would 
decay into photon pairs with a narrow width. From another side, studies of prompt dipho- 
ton production provide a particularly clean test of perturbative Quantum Chromodynamics 
(pQCD) and soft-gluon resummation methods implemented in theoretical calculations since 
the produced photons are largely insensitive to the effects of final-state hadronization, and 
their energies and directions can be measured with high precision. Therefore, it is essential to 
have an accurate pQCD predictions for corresponding cross sections and related kinematical 
distributions. 

The leading contributions to the production of prompt photon pairs at hadron collisions 
are the quark-antiquark annihilation qq — > 77, gluon-gluon fusion gg — > 77 and quark-gluon 
scattering qg — > 775 subprocesses. Prompt photons may also originate from single or double 
fragmentation processes of the partons produced in the hard interaction [11]. However, a 
photon isolation requirement which involved in the measurements [IHS] significantly reduces 
the rate for these processes^. In the framework of standard QCD, theoretical calculations 
of diphoton production cross sections have been carried out at next-to-leading order (NLO) 
[T2] and next-to-next-to-leading-log approximation (NNLL) [13]. A fixed-order NLO pQCD 
calculations [12] implemented by the diphox program accounts fragmentation subprocesses, 
but the gluon-gluon fusion gg — > 77 contribution is considered only at LO. In the NNLL 
calculations [13] implemented by the RESBOS program the effects of initial state soft gluon 
radiation in the NLO calculations are analytically resummed to all orders in the strong 
coupling constant. Both these calculations reproduce the main features of recent Tevatron 
and LHC data [THS], but none of them describes all aspects of the data. So, there is significant 
underestimation [4J of diphoton cross sections measured by the CMS collaboration in the 
regions of phase space where two photons have an azimuthal angle difference A0 77 < 2.8. 
The similar observation was made [5] by the ATLAS collaboration: more photon pairs 
are seen in data at low A</> 77 values, while the NLO and NNLL pQCD predictions favour 
a larger back-to-back production (A0 77 ~ 71"). At the Tevatron energies, there is the same 
problem [IHS] common to both NLO [12] and NNLL [13] pQCD calculations in the description 
of events with low diphoton mass M 77 or low azimuthal angle distance A0 77 (and also with 
moderate diphoton transverse momentum p^). Such disagreement between the collinear 
QCD predictions and the available data indicates (THS] the necessity of including higher 

1 Usually the photons are called prompt if they are coupled to the interacting quarks. 
2 See discussion in Section 2. 
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order corrections beyond NLO. 

We note, however, that an alternative description of diphoton production at high en- 
ergies can be achieved within the framework of /^-factorization QCD approach [HI [15]. 
This approach is based on famous Balitsky-Fadin-Kuraev-Lipatov (BFKL) [16] or Ciafaloni- 
Catani-Fiorani-Marchesini (CCFM) [17] equations and provides solid theoretical grounds for 
the effects of initial state gluon radiation and intrinsic non-zero parton transverse momen- 
tum. A detailed description and discussion of the /c^-factorization formalism can be found, 
for example, in reviews [18J. Here we would like to only mention that the main part of high- 
order radiative QCD corrections is naturally included into the leading-order /c^-factorization 
formalism. Moreover, the soft gluon resummation formulas implemented in NNLL calcula- 
tions [13] are the result of the approximate treatment of the solutions of CCFM equation, as 
it was shown in [19]. The /^-factorization approach has been successfully applied recently, in 
particular, to describe the prompt photon and Drell-Yan pair production at HERA, Tevatron 
and LHC [2UH25]. 

In the present paper we apply the /cT-factorization QCD approach to diphoton production 
in pp and pp collisions at high energies. Our main goal is to give a systematic analysis of 
available Tevatron data taken by the D0 and CDF collaborations [IH3] and recent LHC 
measurements [HE] performed by the CMS and ATLAS collaborations. The consideration is 
based on the off-shell (depending on the non-zero transverse momenta kx of incoming quarks 
and gluons) production amplitudes of quark-antiquark annihilation qq — > 77, quark-gluon 
scattering qg — >■ 775 and gluon-gluon fusion gg — > 77 subprocessefl The unintegrated (* T - 
dependent) quark and gluon densities in a proton are defined by using the Kimber-Martin- 
Ryskin (KMR) prescription [26]. This approach is the simple formalism to construct the 
unintegrated parton distributions from the known conventional ones. We calculate total and 
differential diphoton production cross sections and estimate the theoretical uncertainties of 
our predictions. Special attention is put to the specific kinematic properties of the produced 
photon pair which are strongly related to the non-zero transverse momenta of initial partons. 
Such investigations in the framework of //^-factorization approach are performed for the first 
time. 

The outline of our paper is following. In Section 2 we recall shortly the basic formulas 
of the fcr-factorization approach with a brief review of calculation steps. In Section 3 we 
present the numerical results of our calculations and a discussion. Section 4 contains our 
conclusions. 

2 Theoretical framework 

As it was mentioned above, the leading contributions to the diphoton production at high 
energies are the quark-antiquark annihilation, quark-gluon scattering and gluon-gluon fusion 

3 In the case of gluon-gluon fusion we will neglect the transverse momenta of incoming gluons in the 
corresponding production amplitude but keep true off-shell subprocess kinematics. See Section 2 for more 
details. 
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subprocesses: 

q(h) + q(k 2 ) ^7(^3) + 7(^4), (1) 

q(h) + g(k 2 ) ->• 7(^3) + 7(^4) + g(fc 5 ), (2) 
g(k 1 )+g(k 2 )->'y(k 3 )+'y(k 4 ), (3) 

where the four-momenta of all particles are given in parentheses. Although last two of them 
are higher-order processes, their contributions are quantitatively comparable to those from 
quark- ant iquark annihilation in the diphoton invariant mass range of interest, due to the 
significant gluon luminosity in this kinematical region. 

Let us start from the kinematics. In the center-of-mass frame of colliding particles we 
can write: 



P (1) = ^(1,0,0,1), p^ = ^(1,0,0,-1), (4) 

where p^ and p < - 2 ' ) are the four-momenta of colliding protons, y/s is the total energy of 
the process under consideration and we neglect the masses of protons. The initial parton 
four-momenta in the high energy limit can be written as follows: 

k\ = xip^ + k 1T , k 2 = x 2 p^ + k 2 T, (5) 

where k\T and k 2 T are their transverse four-momenta. It is important that kf T = — k\ T 7^ 0, 
k|r = ~k 2 T 7^ 0. From the conservation laws we can easily obtain the following relations for 
subprocesses (1) and (3): 

kiT + k 2T = k 3T + k 4T , (6) 
x^s = \k 3T \e y3 + \k 4T \e y \ (7) 
x 2V ^ = |k 3T |e-^ + \k 4T \e- y \ (8) 
and the similar ones for subprocess (2): 

ki T + k 2T = k 3T + k 4T + k 5T , (9) 

x^s = \k 3T \e yz + |k 4T |e^ 4 + m bT e y '% (10) 

x 2 Vs = \k 3T \e~ y3 + |k 4T |e- y4 + m hT e~ y \ (11) 

where k 3 T, ^4T and k§T are the transverse four-momenta of corresponding particles, y 3 , y 4 
and 1/5 are their center-of-mass rapidities and m^r is the transverse mass of produced quark 
(or antiquark) having mass m. 

The off-shell partonic amplitudes of (1) and (2) can be written as follows: 



xu(k 2 



M(q*q* ->■ 77) = e 2 e 2 e^k 3 )e u (k 4 ) x 

„,„ ki-k + m h-h + m 

(ki — k 3 ) 2 — m 2 (ki — k^) 2 — m 2 



u(ki), 



(12) 



M(q*g* ^ liq ) = e 2 e 2 q gt a e x {k 2 )e,{h)e u {k,)Y,n V \ (13) 



n 



J~ 9 



■fii/X 



u(hh l 
u(fa)Y 
u(fa)Y 



fa + k$ + m 
(fa + k 5 ) 2 — m 2 



fa + k 2 + m 
(ki + k 2 ) 2 — m A 



k 5 + m 



fa^ 



r7 



(fa T K 5 

fa + fa + m 
{k 4 + fa) 2 — m 2 1 (ki — fa) 2 — m' 



r7 



fa + fa + m 
(fa + fa) 2 — 

fa — fa + m 



puX 



J~ 5 



-pfivX 



u(fah x '"' 

u(fa)-f 



k 2 + m 



(fa - fa) 2 - m 
fa + fa + m 



(fa + fa 



u(fa)Y 



fa 



- 77T 

k 2 + m 



r7 



fa — fa + m 
(fa — ki) 2 — m 2 

fa — ki + m 



r7" 



(fa- 
fa 



fa) 2 



r7M*i). 
7"ti(A;i), 

7*«(*i). 



(14) 
(15) 
(16) 
(17) 



fc 3 + m 



fcO 2 



-,^u(k 



(19) 



(fa-fa) 2 -m? (Ali-res; 

where e and e g are the electron and quark (fractional) electric charges, g is the strong charge, 
e(k 2 ), e(fa) and £(£4) are the polarization four- vectors of corresponding particles and a is 
the eight-fold color index. When we calculate the matrix elements squared, the summation 
on the polarizations of produced photons is carried out by usual covariant formula: 



flU 



^(h)e^(h) = -g 
where i = 3,4. In the case of initial off-shell gluon we apply the BFKL prescription [1 



Y,nfaY v (fa) 



h " 2 T hj 2T 

~k 2 ' 

2T 



(20) 



(21) 



This formula converges to the usual one (20) after azimuthal angle averaging in the |k 2 <r| — > 
limit. Since we do not neglect the transverse momentum of incoming quarks, the standard 
on-shell quark spin density matrix has to be replaced by a more complicated expression. 
To evaluate it we will follow an approximation proposed in [20J. We "extend" the original 
diagram and consider the off-shell quark line as internal line in the extended diagram. The 
extended process looks like follows: the initial on-shell quark with four-momentum p and 
mass m radiates a quantum (say, photon or gluon) and becomes an off-shell quark with 
four-momentum k. So, for the extended diagram squared we can write: 



\M\ 



tr 



7 M — -7 u(p)u(p)^ v — -Ta 



k 2 — m 2 



k 2 — m 2 M 



(22) 



where T is the rest of the original matrix element. The expression presented between T 7 ' 
and Tjj, now plays the role of the off-shell quark spin density matrix. Using the standard 
on-shell condition u(p)u(p) = p + m and performing the Dirac algebra one obtains in the 
massless limit m — > 0: 



\M\ 



(k 



2\2 



tr 



T»(k 2 p-2(p-k)k)% 



(23) 
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Now, using the Sudakov decomposition k = xp + (where kx is the off-shell quark non- 
zero transverse four-momentum, k 2 = k\ = — k\) and neglecting the second term in the 
parentheses in (23) in the small- x limit, we easily obtain: 

\M\ 2 ~^Lt x [T»xp%]. (24) 

Essentially, we have neglected here the negative light-cone momentum fraction of the incom- 
ing quark. The properly normalized off-shell spin density matrix is given by xp, while the 
factor 2/xk 2 has to be attributed to the quark distribution function (determining its leading 
behavior). With this normalization, we successfully recover the on-shell collinear limit when 
k is collinear with p. Further calculations are straighforward and in other respects follow 
the standard QCD Feynman rules. We only mention that the method of orthogonal ampli- 
tudes [27] has been applied to avoid the long output. The evaluation of traces in (12) - 
(19) was done using the algebraic manipulation system FORM [25] . 

The matrix element squared of gluon-gluon fusion subprocess (3) was calculated a long 
time ago in the on-shell limit |kiy| — > 0, |k2r| — > 0. The simple analytical expression can 
be found, for example, in [29]. In our phenomenological study, we apply it with one remark: 
numerically, we keep the exact off-shell kinematics given by (5) — (8). 

In according to the /^-factorization theorem, to calculate the cross section of diphoton 
production one should convolute off-shell partonic cross sections with the relevant uninte- 
grated quark and/or gluon distributions in a proton: 



<?= ^2 J £ii(si>E 2 ,k? r ,]4r)/i(ai,k^ 



(25) 



h3=q>9 



where a*j(xi,X2, kf T , k 2 T ) is the off-shell partonic cross section and fi(x, k^, /i 2 ) is the unin- 
tegrated parton densities in a proton. The contributions from quark- ant iquark annihilation, 
quark-gluon scattering and gluon-gluon fusion can be easily rewritten as follows: 



xf q (x u k 2 T , /J 2 )f q (x 2 , k 2 T , /j 2 )dk 2 1T dk 2 2T dkl T dy 3 dy 4 



(26) 



2tt 2tt 



a 



^/ 2 56,3(U) 2 |A1(9V " 779)|2X 



(27) 
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2vr 2vr 2vr 2vr 



a= / Tr~T~ dM^^77)| 2 x 

J 167r(a;ia;2S i )^ 

xf g (x 1: k\ T , n 2 )f g (x 2 , kj T , ii 2 )dk\ T dkl T dkl T dy z dy±- 



(2f 



2tt 2tt ' 

where <pi, <f) 2 , 4>3 and 04 are the azimuthal angles of initial partons and produced photons, 
respectively. If we average these expressions over <pi and <j) 2 and take the limit |kiy| — > and 
\k 2 j-\ — > 0, then we recover the corresponding formulas in the collinear QCD factorization. 
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Concerning the unintegrated quark and gluon densities in a proton, we apply the Kimber- 
Martin-Ryskin (KMR) approach to calculate them. The KMR approach is the formalism 
to construct the unintegrated parton distributions from the known conventional ones. In this 
approximation the unintegrated quark and gluon distributions are given by 



r 2 ,,2\ _ ; ,1,2 ,.2\ a s(^r) 



fq( x i ^Ti A* ) ~~ ^o(ky, /i 



2tt 



X 



1 

x J dz 

x 



(29) 



P„(z)-q (-Xt) Q(A-z) + P gg (z)-g (-,k 2 

Z \ Z / z \ z 



fg(X,]i T ,fi ) —Tg(k T ,fl ^ 



2n 



± 

[ dz ^P„(z)^q(^)+P„(z)^g(^)e(A-z) 

J Z V Z ' Z V Z ' 



(30) 



where P a b(z) are the usual unregulated LO DGLAP splitting functions. The theta functions 
which appears in (29) and (30) imply the angular-ordering constraint A = fi/(fJL+ \kr\) specif- 
ically to the last evolution step to regulate the soft gluon singularities. Numerically, for the 
input we have used leading-order parton densities xq(x,fi 2 ) and xg(x,fi 2 ) from MSTW'2008 
set [30]. 

We note that perturbation theory becomes nonapplicable when the wavelength of the 
emitted photon (in the emitting quark rest frame) becomes larger that the typical hadronic 
scale 0(1 GeV -1 ). Then the nonperturbative effects of hadronization or fragmentation must 
be taken into account. Accordingly, the calculated cross section can be split into two pieses 



O = Odir(Ai 2 ) + CT fragm (/X 



(31) 



with o"dir(/i 2 ) representing the perturbative contribution and cr fragm (/i 2 ) the fragmentation 
contribution. In the fragmentation processes photons are produced through the fragmen- 
tation of a parton (produced in a hard subprocess) into a single photon carrying a large 
fraction z of parent parton momentum. These processes are described in terms of quark-to- 
photon D q -> 7 (z, /i 2 ) and gluon-to-photon D g ^(z, /i 2 ) fragmentation functions [31 J . In our 
calculations for quark-gluon scattering (2) we choose the fragmentation scale fi 2 to be the 
invariant mass of the quark + photon subsystem, and restrict o"dir(yU 2 ) to \i > M ~ 1 GeV. 
Under this condition, the contribution <y^ v {^ 2 ) is free from divergences (so the mass of the 
light quark m can be safely sent to zero) and we checked that the sensitivity of our results 
to the choice of M is reasonably soft. As far as the fragmentation contribution crf ragm (/z 2 ) is 
concerned, its size is dramatically reduced by the photon isolation cuts implemented in the 
experimental analysis [THS]. The isolation condition required that the hadronic transverse 
energy E^ ad , deposited inside a cone with aperture R centered around the photon direction 
in the pseudo-rapidity rj and azimuthal angle plane, is smaller than some value E™ x - In 
the recent measurements [IHH] performed at the Tevatron and LHC, these parameters were 
taken R ~ 0.4 and i?™ ax ~ 1 — 2 GeV. According to the estimates [3], the contribution from 
Cfragm(A* 2 ) amounts to about 10 — 15% of the visible cross section and therefore is neglected 
in our consideration. 
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The multidimensional integrations in (26) — (28) have been performed by the means of 
Monte Carlo technique, using the routine VEGAS |32j. The full C++ code is available from 
the author on request. 

3 Numerical results 

We now are in a position to present our numerical results. First we describe our theoret- 
ical input and the kinematic conditions. After we fixed the unintegrated gluon distributions, 
the cross sections (26) — (28) depend on the renormalization and factorization scales lir 
and ix p. Numerically, we set them to be equal to n R = fx F = (£M 77 ) 2 . In order to estimate 
the scale uncertainties of our calculations we vary the parameter £ between 1/2 and 2 about 
the default value £ = 1. Since the expression for the off-shell quark spin density matrix has 
been derived in the massless approximation, numerically we neglect the quark masses. We 
use the LO formula for the strong coupling constant a s (/i 2 ) with rif = 4 active quark flavors 
at Aqcd = 200 MeV, so that a s (M|) = 0.1232. 

Experimental data for prompt photon pair production in pp collisions at the Tevatron 
come from both the CDF [U[2] and D0 |3J collaborations. Several differential cross section 
have been determined: as a function of the diphoton invariant mass M 77 , the transverse 
momentum and rapidity of photon pair p^ and y 77 , the azimuthal angle difference between 
the produced photons A0 77 and the cosine of the polar scattering angle 9* of leading photon 
in the Collins-Soper frame. Recent D0 data [3] refer to the kinematic region defined by |r7 71 1 < 
0.9, |rf 2 | < 0.9, pf > 21 GeV and pf > 20 GeV with the total energy yfe = 1960 GeV, 
and an additional cut M 77 > pjJ was applied. The measurements of double differential cross 
sections da / dM^dpTp , da / 'dM 77 dA0 77 and da / dM 11 d\ cos#*| for three subdivisions of M 77 
range (namely 30 < M 77 < 50 GeV, 50 < M 77 < 80 GeV and 80 < M 77 < 350 GeV) have 
been presented also. The CDF data [2] refer to the kinematic region defined by \y ll \ < 1, 
\y l2 \ < 1, p T l > 17 GeV and p^ > 15 GeV. In this analysis, the diphoton cross section as 
a function of variable z = p 1 ^ jp 1 ^ , the ratio of sub- leading to leading photon transverse 
momentum (0 < z < 1), has been measured additionally. Also two specific kinematic cases 
have been examined separately: differential cross sections for M 77 > pip and M 77 < p7p ' . 
Previous CDF data pL| refer to the kinematic region defined by |r/ 7l | < 0.9, |rf 2 | < 0.9, 
pip > 14 GeV, pf > 13 GeV and y/s = 1960 GeV. 

First experimental data for prompt diphoton production in pp collisions at the LHC come 
from the CMS jl] and ATLAS |5J collaborations. The CMS data [I] have been obtaned in 
the kinematic region defined by |r/ 7l | < 2.5, |r/ 72 | < 2.5 (except pseudo-rapidity region 
1.44 < 77 < 1.57 for both photons), p~fi > 23 GeV and pf > 20 GeV with the total 
energy y/s = 7 TeV. The measurements for more tight central region with |r] 7l | < 1.44 and 
|t7 72 | < 1.44 have been presented additionally. The ATLAS data j5] refers to the kinematic 
region defined by |?7 7l | < 2.37, |?7 72 | < 2.37 (with the exclusion of the region 1.37 < 77 < 1.52 
for both photons), p T x > 16 GeV and pl^ > 16 GeV with the with the same total energy. 

The results of our calculations are shown in Figs. 1 — 10 in comparison with the Tevatron 
data [IH3] and in Figs. 11 — 15 in comparison with the LHC data [HE]. The solid histograms 
are obtained by fixing both the factorization and renormalization scales at the default value, 
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Source 


a(pp -»• 77 + X) [pb] 


/^-factorization (KMR) 
NLO pQCD [12] (diphox) 
NNLL pQCD [13J (resbos) 
CDF data [2] 


8.87±5:|| (scales) 
10.58 ±0.55 
11.31 ±2.45 
12.47 ±0.21 (stat.) ±3.74 (syst.) 



Table 1: The total cross section of diphoton production in pp collisions calculated in the 
kinematical region defined by \y^\ < 1, \y~< 2 \ < 1, p^ 1 > 17 GeV and pf > 15 GeV at 
y/s = 1960 GeV. 

whereas the upper and lower dashed histograms correspond to the scale variation as it was 
described above (in left panels). Also, in Figs. 1 — 6 we plot for comparison the NNLL pQCD 
predictions [13J (as given by resbos program) listed in [HE]. The relative contributions to 
the diphoton cross sections from quark-antiquark annihilation, quark-gluon scattering and 
gluon-gluon fusion subprocesses are shown separately by the dash-dotted, dotted and dashed 
histograms, respectively, in right panels. 

We note that different kinematic variables under consideration probe different aspects of 
the diphoton production mechanism. For instance, the M 77 spectrum is particularly sensitive 
to potential contributions from physics beyond the SM. The cos 9* distribution probes the 
angular momentum of the final state, which should be different for QCD-mediated production 
as compared, for example, to the decay of a scalar Higgs boson H — > 77 (see [T3]). 

We find that the /^-factorization approach reasonably well describe a full set of experi- 
mental data taken at the Tevatron and LHC. Moreover, the shape and absolute normalization 
of measured cross sections are adequately reproduced within the theoretical and experimen- 
tal uncertainties. At the Tevatron, there is some discrepancy between recent CDF data [2] 
and our predictions at very high diphoton invariant masses (see Fig. 7). Note, however, 
that the small-x approximation used in our consideration is broken at high M 77 and the 
more accurate treatment of, in particular, off-shell quark spin density matrix is needed. 
This point, of course, is of importance but it is out of our present consideration. From 
another side, we would like to point out that /cT-factorization approach provides a better 
description of all kinematic variables (in comparison with the NLO/NNLL pQCD calcula- 
tions p2J[l3]) at low M 77 values, where expected effects connected with the small-a: physics 
becomes important. In the collinear QCD factorization, observed large discrepancy between 
data and NLO/NNLL pQCD predictions p^2lfT3] at low M 77 attributes to a fragmentation 
contributions and higher-order QCD corrections to the gluon-gluon fusion subprocess [TH5]. 
However, the main part of these corrections are already effectively taken into account in 
our calculations. It is a radiative QCD corrections described by the ladder-type diagrams 
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Source 


a(pp ->• 77 + X) [pb] 


/cT-factorization (KMR) 
NLO pQCD pj] (diphox) 
CMS data 


29.9i£;2 (scales) 
27.3 ±i;g (scales) ±1.1 (PDFs) 
31.0 ±1.8 (stat.) ^-o ( syst _) ±12 (lumi.) 



Table 2: The total cross section of diphoton production in pp collisions calculated in the 
kinematical region defined by |^| < 1.44, |rf 2 | < 1.44, p^ > 23 GeV and pf > 20 GeV at 
y/s = 7 TeV. 

which are naturally included into the fc^-factorization formalism at LO levefl. At the LHC, 
our predictions agree well with the CMS data [4] in a whole kinematical range and slightly 
overestimate the ATLAS data [5] at low values of A0 77 and high values of p^f . This point, 
probably, can indicate the inconsistecy of the CMS and ATLAS data and, of course, needs 
in a further theoretical and experimental investigations. Note also that the NLO pQCD 
predictions [12] underestimate the measured cos 9* distributions at both the Tevatron and 
LHC energies [2H5]. This underestimation is more significant for the central rapidity range 
(\r]\ < 1.44) at the LHC. Our predictions agree well with the data [2H5]. It could be im- 
portant for experimental detection of Higgs boson signal and further investigations of Higgs 
properties. 

As it was expected, the main important properties of the /c^-factorization approach clearly 
manifest themselves in the calculated pj 7 , z and A0 77 distributions. In the naive LO pQCD 
approximation, photons produced in 2 — > 2 subprocesses are back-to-back in the trans- 
verse plane and are balanced in pt due to momentum conservation. Therefore, these three 
ditsributions must be simply a delta functions: da/dplp ~ 8(p^P), da/dz ~ 8(z — 1) and 
da I g?A0 77 ~ <5(A0 77 — 7r). When higher-order QCD processes are considered, the pres- 
ence of additional quarks and/or gluons in the final state allows these distributions to be 
more spread. In the /^-factorization approach, taking into account the non- vanishing initial 
parton transverse momentum leads to the violation of back-to-back kinematics even at LO 
approximation. So, despite the fact that one from contributed subprocesses is the 2—7-3 
subprocess (that makes the difference between the A>r-factorization predictions and ones from 
the collinear approximation of QCD not well pronounced), obtained perfect description of 
measured A0 77 distributions |TH5| is notable. Specially we point out a reasonable description 
of the Tevatron and LHC data at low A0 77 values where the significant underestimation of 
the data by the NLO/NNLL pQCD calculations p2J[T3] was observed (see also Figs. 2, 4 
and 6). The important role of such angular correlations for understanding an interaction 
dynamics is well known [18] . In particular, as it was demonstrated in [33] , these correlations 

4 See review [TB] for more details. 
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Source 


a(pp -»• 77 + X) [pb] 


/cT-factorization (KMR) 
NLO pQCD pj] (diphox) 
CMS data 


59.2™ (scales) 
52.7 ±£;1 (scales) ±2.0 (PDFs) 
62.4 ± 3.6 (stat.) t\ j (syst.) ±2.5 (lumi.) 



Table 3: The total cross section of diphoton production in pp collisions calculated in the 
kinematical region defined by |?7 7l | < 2.5, |?7 72 | < 2.5 (except pseudo-rapidity region 1.44 < 
T] < 1.57 for both photons), p£ > 23 GeV and pf > 20 GeV at y/s = 7 TeV. 

strongly depend on the unintegrated parton densities involved in the calculations and can 
be used as a test to distinguish the different approaches of non-collinear parton evolution^. 
We would like to note also that the scale uncertainties of our calculations are quite small at 
A0 77 ~ and increases when A0 77 — >■ it. Similar to A0 77 distributions, the calculated p7p 
distributions also directly connected with the unintegrated quark and gluon distributions due 
to transverse momentum conservation. Perfect overall agreement of our calculations and the 
data on the A0 77 and p'P shows that the KMR prescription for evaluation of unintegrated 
parton densities reproduces well the transverse momenta of initial quarks and gluons in a 
probed kinematical region. 

Special interest for the fc^-factorization phenomenology is connected also with the events 
with low invariant mass and high transverse momentum of produced photon pair. The CDF 
collaboration has studied this kinematical region using a special condition M 77 < pf in their 
experimental analysis [2]. Within the mentioned above cuts on the transverse momenta of 
the produced photons, such measurements provide us with an additional test for unintegrated 
quark and gluon distributions f a ( ) at moderate and high k|, values. The events with 

back-to-back kinematics are strongly suppressed in this case (see Fig. 10). In contrast with 
the NLO/NNLL pQCD calculations [12j[T3], the /cy-factorization predictions agree well with 
the CDF data. An opposite condition (namely, M 77 > pjJ) corresponds to the selection of 
events with kinematics similar to the decay of a heavy particle with low transverse momentum 
into a photon pair, such as, for example, H — > 77. Such events are better described [2] by 
the NLO pQCD calculations [12] than ones with low invariant mass and high transverse 
momentum. Our predictions agree within the uncertainties with the data [2], but overall 
agreement is a bit worse. 

Concerning a relative contributions to the calculated diphoton production cross sections, 
we find that quark- ant iquark subprocess dominates at both the Tevatron and LHC energies 

5 We note, however, that at moment the KMR prescription is only one provide us with the unintegrated 
quark densities f q {x, , /1 2 ) which can be used in a wide region of x and k^ in phenomenological studies. 
Thus, the problem of estimation of theoretical uncertainties connected with unintegrated parton densities is 
still open. 
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in the probed kinematical regional. The contribution from quark-gluon scattering is similar 
to gluon-gluon fusion one at low diphoton invariant masses and overshoot gluon-gluon fusion 
contribution at high M 77 values. As it was expected [HE], the role of both these subprocesses 
increases significantly at the LHC energies. We note also that the contribution from qq — > 
77$ subprocess, which is order of 0(a s a 2 ), is already taken into account in our calculations. 
To be precise, this contribution is effectively included to the quark-antiquark annihilation one 
due to the initial state gluon radiation (see also [IE])- Of course, it is in a clear contrast with 
the collinear QCD factorization where these contributions should be calculated separately. 

The predicted total cross sections for the Tevatron conditions are listed in Table 1 and 
for LHC conditions are listed in Tables 2 and 3, respectively. Our results are slightly below 
the Tevatron data and NLO/NNLL pQCD predictions [T2l[T3] but agree with them within 
the theoretical and experimental uncertainties. However, at the LHC energy, our predictions 
overshoot the NLO pQCD ones [12] and are more close to the data. It could be connected 
with the small- a; effects which becomes more important for the LHC conditions. 

In conclusion, we would like to stress a number of important achievements demonstrated 
by the /c^-factorization QCD approach. As a general feature, the /c^-factorization predictions 
are found to be perfectly compatible with the available data, in particular, on the production 
of prompt photons, Drell-Yan lepton pairs, heavy quarks as well as various quarkonium states 
at modern colliders [2014251 [33] [35 | 136] . The fc T -factorization approach succeeds in describing 
the polarization phenomena observed in ep and pp interactions (see, for example, [35*36] and 
references therein). The latter is essentially related to the initial gluon off-shellness, which 
dominates the gluon polarization properties and has a considerable impact on the kinematics. 
So, we believe that the /^-factorization formalism holds a possible key to understanding the 
production dynamics at high energies. 

4 Conclusions 

We have investigated the isolated prompt photon pair production in pp and pp collisions at 
the Tevatron and LHC energies within the framework of the /^-factorization QCD approach. 
Our consideration is based on the quark-antiquark annihilation, quark-gluon scattering and 
gluon-gluon fusion subprocesses, where the non-zero transverse momenta of incoming partons 
are taken into account. The unintegrated parton densities in a proton are determined using 
the Kimber-Martin-Ryskin prescription. We obtained a reasonable well agreement between 
our predictions and recent data taken by the D0, CDF, CMS and ATLAS collaborations. The 
achieved description of the data in several kinematical regions (in particular, at low azimuthal 
angle distance between produced photons) is better as compared to the one obtained in 
the framework of collinear QCD factorization (with NLO or even NNLL accuracy). Also, 
we describe reasonably well the Tevatron and LHC data on the cos 8* distributions of the 
diphoton production. Such quantities (underestimated by the NLO pQCD calculations) 
probe the angular momentum of the final state and different for QCD-mediated production 
as compared to the decay of a scalar Higgs boson. It could be important for experimental 
detection of Higgs and further theoretical and experimental investigations in this field. We 

6 We strongly disagree with the estimations [34] of relative contributions from the gluon-gluon fusion and 
quark-antiquark annihilation subprocesses performed for the Tevatron energy. 
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find that the scale uncertainties of our calculations are quite small at A0 77 ~ 0, although 
they are exceed, of course, the uncertainties of NLO/NNLL pQCD calculations in general. 
Moreover, the problem of estimation of our theoretical uncertainties (connected in part with 
the non-collinear parton evolution scheme) is still open and further theoretical attempts (in 
order to investigate the unintegrated quark distributions more detail) are necessary to reduce 
this uncertainty. It is important for further studies of small-rr physics at hadron colliders, 
and, in particular, for searches of effects of new physics beyond the SM at the LHC. 
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Figure 1: The differential cross section of prompt photon pair production in pp collisions 
at the Tevatron as a function of diphoton invariant mass M 77 and transverse momentum 
p7p . Left panel: the solid histograms correspond to the results obtained with the KMR 
parton densities at the default scale, and the upper and lower dashed histograms correspond 
to standard scale variations, as it is described in the text. The dotted histograms repre- 
sent the NNLL pQCD predictions [T3] (as given by RESBOS program) listed in [3J. Right 
panel: the different contributions to the diphoton cross sections calculated at the default 
scale. The dashed, dash-dotted and dotted histograms correspond to the gluon-gluon fusion, 
quark- ant iquark annihilation and quark-gluon scattering subprocesses. The solid histograms 
represent the sum of these contributions. The experimental data are from D0 [3j. 



16 




A<t>w[rad] 



S3 

a. 



CD 
CO 

O 1 

_o 10 1 
■a 

D 
"O 



— i — i — i — i — i — i — i — 




10' 



0.0 



—I I I I I I 1_ 



DO ► 



0.2 



0.4 



0.6 

|cos e*| 



T3 
nj 

D. 10 1 

1 10° 



10" 



10"' 



T 1 1 ' 


I 1 1 1 1 

DO i-"- 




1 , '■ 




: ; \ 




-■ 


■ 




J I I 





A^[rad] 




Figure 2: The differential cross section of prompt photon pair production in pp collisions at 
the Tevatron as a function of A0 77 and | cos#*|. Notation of all histograms is the same as 
in Fig. 2. The experimental data are from D0 [3]. 



17 



10' 



30 



50 



P7 [GeV] 



IU 



60 80 

P7 [GeV] 



1 — ■ 


— 1 1 


i — i — i — i — i — i — 

DO 

80 < M yy < 350 GeV 


1 1 1 









20 40 60 80 100 



Figure 3: The double differential cross section da / dM yy dp^f of prompt photon pair produc- 
tion in pp collisions at the Tevatron calculated for three subdivisions of M 77 range. The solid 
histograms correspond to the results obtained with the KMR parton densities at the default 
scale, and the upper and lower dashed histograms correspond to standard scale variations, as 
it is described in the text. The dotted histograms represent the NNLL pQCD predictions [T3] 
(as given by RESBOS program) listed in [3]. The experimental data are from D0 [3]. 
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Figure 4: The double differential cross section da / dM 71 dA(f)^ of prompt photon pair produc- 
tion in pp collisions at the Tevatron calculated for three subdivisions of M 77 range. Notation 
of all histograms is the same as in Fig. 3. The experimental data are from D0 |3J. 
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production in pp collisions at the Tevatron calculated for three subdivisions of M 77 range. 
Notation of all histograms is the same as in Fig. 3. The experimental data are from D0 [3J. 



20 




10' 



A^[rad] 



10"' 



A^[rad] 



Figure 6: The differential cross section of prompt photon pair production in pp collisions at 
the Tevatron. Notation of all histograms is the same as in Fig. 2. The experimental data 
are from CDF [TJ. 
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Figure 7: The differential cross section of prompt photon pair production in pp collisions at 
the Tevatron as a function of diphoton invariant mass M 77 . Two kinematic cases are shown 
separately: differential cross sections for M 77 > pf and M 77 < pip . The solid histograms 
correspond to the results obtained with the KMR parton densities at the default scale, and 
the upper and lower dashed histograms correspond to standard scale variations, as it is 
described in the text. The experimental data are from CDF [2]. 
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Figure 8: The differential cross section of prompt photon pair production in pp collisions at 
the Tevatron as a function of diphoton transverse momentum p^J. Two kinematic cases are 
shown separately: differential cross sections for M 77 > p7p and M 77 < p7p . Notation of all 
histograms is the same as in Fig. 7. The experimental data are from CDF [2]. 
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Figure 9: The differential cross section of prompt photon pair production in pp collisions at 
the Tevatron as a function of diphoton rapidity y 77 and variable z = pj? /p^ ■ Two kinematic 
cases are shown separately: differential cross sections for M 77 > pl^ and M 77 < pip . Notation 
of all histograms is the same as in Fig. 7. The experimental data are from CDF [2]. 
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Figure 10: The differential cross section of prompt photon pair production in pp collisions at 
the Tevatron as a function of A0 77 and cos#*. Two kinematic cases are shown separately: 
differential cross sections for M 77 > pp and M 77 < pip . Notation of all histograms is the 
same as in Fig. 7. The experimental data are from CDF [2]. 
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Figure 11: The differential cross section of prompt photon pair production in pp collisions 
at the LHC as a function of diphoton invariant mass. Left panel: the solid histograms cor- 
respond to the results obtained with the KMR parton densities at the default scale, and the 
upper and lower dashed histograms correspond to standard scale variations, as it is described 
in the text. Right panel: the different contributions to the diphoton cross sections calcu- 
lated at the default scale. The dashed, dash-dotted and dotted histograms correspond to the 
gluon-gluon fusion, quark-antiquark annihilation and quark-gluon scattering subprocesses. 
The solid histograms represent the sum of these contributions. The experimental data are 
from CMS 0]. 
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Figure 12: The differential cross section of prompt photon pair production in pp collisions 
at the LHC as a function of diphoton transverse momentum. Notation of all histograms is 
the same as in Fig. 11. The experimental data are from CMS [I]. 
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Figure 13: The differential cross section of prompt photon pair production in pp collisions 
at the LHC as a function of A0 77 . Notation of all histograms is the same as in Fig. 11. The 
experimental data are from CMS [4J. 
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Figure 14: The differential cross section of prompt photon pair production in pp collisions 
at the LHC as a function of | cos8*\. Notation of all histograms is the same as in Fig. 11. 
The experimental data are from CMS [I]. 
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Figure 15: The differential cross section of prompt photon pair production in pp collisions 
at the LHC Notation of all histograms is the same as in Fig. 11. The experimental data are 
from ATLAS 0. 
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